Impact of mutations on the stability of SARS-CoV-2 nucleocapsid protein structure.

The nucleocapsid (N) protein of SARS-CoV-2 is known to participate in various host cellular processes, including interferon inhibition, RNA interference, apoptosis, and regulation of virus life cycles. Additionally, it has potential as a diagnostic antigen and/or immunogen. Our research focuses on examining structural changes caused by mutations in the N protein. We have modeled the complete tertiary structure of native and mutated forms of the N protein using Alphafold2. Notably, the N protein contains 3 disordered regions. The focus was on investigating the impact of mutations on the stability of the protein's dimeric structure based on binding free energy calculations (MM-PB/GB-SA) and RMSD fluctuations after MD simulations. The results demonstrated that 28 mutations out of 37 selected mutations analyzed, compared with wild-type N protein, resulted in a stable dimeric structure, while 9 mutations led to destabilization. Our results are important to understand the tertiary structure of the N protein dimer of SARS-CoV-2 and the effect of mutations on it, their behavior in the host cell, as well as for the research of other viruses belonging to the same genus additionally, to anticipate potential strategies for addressing this viral illness․.

Effect of Colchicine Binding Site Inhibitors on the Tubulin Intersubunit Interaction.

Microtubules are dynamic, non-covalent polymers consisting of α- and ß-tubulin subunits that are involved in a wide range of intracellular processes. The polymerization and dynamics of microtubules are regulated by many factors, including small molecules that interact with different sites on the tubulin dimer. Colchicine binding site inhibitors (CBSIs) destabilize microtubules and inhibit tubulin polymerization, leading to cell cycle arrest. Because of their therapeutic potential, the molecular mechanism of CBSI function is an area of active research. Nevertheless, important details of this mechanism have yet to be resolved. In this study, we use atomistic molecular dynamics simulations to show that the binding of CBSIs to the tubulin heterodimer leads to the weakening of tubulin intersubunit interaction. Using atomistic molecular dynamics simulations and binding free energy calculations, we show that CBSIs act as protein-protein interaction inhibitors and destabilize interlinkage between α and ß subunits, which is crucial for longitudinal contacts in the microtubule lattice. Our results offer new insight into the mechanisms of microtubule polymerization inhibition by colchicine and its analogs.

A Study of a Protein-Folding Machine: Transient Rotation of the Polypeptide Backbone Facilitates Rapid Folding of Protein Domains in All-Atom Molecular Dynamics Simulations.

Molecular dynamics simulations of protein folding typically consider the polypeptide chain at equilibrium and in isolation from the cellular components. We argue that in order to understand protein folding as it occurs in vivo, it should be modeled as an active, energy-dependent process, in which the cellular protein-folding machine directly manipulates the polypeptide. We conducted all-atom molecular dynamics simulations of four protein domains, whose folding from the extended state was augmented by the application of rotational force to the C-terminal amino acid, while the movement of the N-terminal amino acid was restrained. We have shown earlier that such a simple manipulation of peptide backbone facilitated the formation of native structures in diverse α-helical peptides. In this study, the simulation protocol was modified, to apply the backbone rotation and movement restriction only for a short time at the start of simulation. This transient application of a mechanical force to the peptide is sufficient to accelerate, by at least an order of magnitude, the folding of four protein domains from different structural classes to their native or native-like conformations. Our in silico experiments show that a compact stable fold may be attained more readily when the motions of the polypeptide are biased by external forces and constraints.

Mechanism of external K+ sensitivity of KCNQ1 channels.

KCNQ1 voltage-gated K+ channels are involved in a wide variety of fundamental physiological processes and exhibit the unique feature of being markedly inhibited by external K+. Despite the potential role of this regulatory mechanism in distinct physiological and pathological processes, its exact underpinnings are not well understood. In this study, using extensive mutagenesis, molecular dynamics simulations, and single-channel recordings, we delineate the molecular mechanism of KCNQ1 modulation by external K+. First, we demonstrate the involvement of the selectivity filter in the external K+ sensitivity of the channel. Then, we show that external K+ binds to the vacant outermost ion coordination site of the selectivity filter inducing a diminution in the unitary conductance of the channel. The larger reduction in the unitary conductance compared to whole-cell currents suggests an additional modulatory effect of external K+ on the channel. Further, we show that the external K+ sensitivity of the heteromeric KCNQ1/KCNE complexes depends on the type of associated KCNE subunits.

Identification of non-classical hCA XII inhibitors using combination of computational approaches for drug design and discovery.

Human carbonic anhydrase XII (hCA XII) isozyme is of high therapeutic value as a pharmacological target and biomarker for different types of cancer. The hCA XII is one of the crucial effectors that regulates extracellular and intracellular pH and affects cancer cell proliferation, invasion, growth and metastasis. Despite the fact that interaction features of hCAs inhibitors with the catalytic site of the enzyme are well described, lack in the selectivity of the traditional hCA inhibitors based on the sulfonamide group or related motifs is an urgent issue. Moreover, drugs containing sulfanomides can cause sulfa allergies. Thus, identification of novel non-classical inhibitors of hCA XII is of high priority and is currently the subject of a vast field of study. This study was devoted to the identification of novel potential hCA XII inhibitors using comprehensive set of computational approaches for drug design discovery: generation and validation of structure- and ligand-based pharmacophore models, molecular docking, re-scoring of virtual screening results with MMGBSA, molecular dynamics simulations, etc. As the results of the study several compounds with alternative to classical inhibitors chemical scaffolds, in particular one of coumarins derivative, have been identified and are of high interest as potential non-classical hCA XII inhibitors.

Combination of consensus and ensemble docking strategies for the discovery of human dihydroorotate dehydrogenase inhibitors.

The inconsistencies in the performance of the virtual screening (VS) process, depending on the used software and structural conformation of the protein, is a challenging issue in the drug design and discovery field. Varying performance, especially in terms of early recognition of the potential hit compounds, negatively affects the whole process and leads to unnecessary waste of the time and resources. Appropriate application of the ensemble docking and consensus-scoring approaches can significantly increase reliability of the VS results. Dihydroorotate dehydrogenase (DHODH) is a key enzyme in the pyrimidine biosynthesis pathway. It is considered as a valuable therapeutic target in cancer, autoimmune and viral diseases. Based on the conducted benchmark study and analysis of the effect of different combinations of the applied methods and approaches, here we suggested a structure-based virtual screening (SBVS) workflow that can be used to increase the reliability of VS.

A new microtubule-stabilizing agent shows potent antiviral effects against African swine fever virus with no cytotoxicity.

African swine fever virus (ASFV) is the causal agent of a fatal disease of domestic swine for which no effective antiviral drugs are available. Recently, it has been shown that microtubule-targeting agents hamper the infection cycle of different viruses. In this study, we conducted in silico screening against the colchicine binding site (CBS) of tubulin and found three new compounds with anti-ASFV activity. The most promising antiviral compound (6b) reduced ASFV replication in a dose-dependent manner (IC50 = 19.5 µM) with no cellular (CC50 > 500 µM) and animal toxicity (up to 100 mg/kg). Results also revealed that compound 6b interfered with ASFV attachment, internalization and egress, with time-of-addition assays, showing that compound 6b has higher antiviral effects when added within 2-8 h post-infection. This compound significantly inhibited viral DNA replication and disrupted viral protein synthesis. Experiments with ASFV-infected porcine macrophages disclosed that antiviral effects of the compound 6b were similar to its effects in Vero cells. Tubulin polymerization assay and confocal microscopy demonstrated that compound 6b promoted tubulin polymerization, acting as a microtubule-stabilizing, rather than a destabilizing agent in cells. In conclusion, this work emphasizes the idea that microtubules can be targets for drug development against ASFV.

Energy-dependent protein folding: modeling how a protein folding machine may work.

Background: Proteins fold robustly and reproducibly in vivo, but many cannot fold in vitro in isolation from cellular components. Despite the remarkable progress that has been achieved by the artificial intelligence approaches in predicting the protein native conformations, the pathways that lead to such conformations, either in vitro or in vivo, remain largely unknown. The slow progress in recapitulating protein folding pathways in silico may be an indication of the fundamental deficiencies in our understanding of folding as it occurs in nature. Here we consider the possibility that protein folding in living cells may not be driven solely by the decrease in Gibbs free energy and propose that protein folding in vivo should be modeled as an active energy-dependent process. The mechanism of action of such a protein folding machine might include direct manipulation of the peptide backbone. Methods: To show the feasibility of a protein folding machine, we conducted molecular dynamics simulations that were augmented by the application of mechanical force to rotate the C-terminal amino acid while simultaneously limiting the N-terminal amino acid movements. Results: Remarkably, the addition of this simple manipulation of peptide backbones to the standard molecular dynamics simulation indeed facilitated the formation of native structures in five diverse alpha-helical peptides. Steric clashes that arise in the peptides due to the forced directional rotation resulted in the behavior of the peptide backbone no longer resembling a freely jointed chain. Conclusions: These simulations show the feasibility of a protein folding machine operating under the conditions when the movements of the polypeptide backbone are restricted by applying external forces and constraints. Further investigation is needed to see whether such an effect may play a role during co-translational protein folding in vivo and how it can be utilized to facilitate folding of proteins in artificial environments.

In silico study of colchicine resistance molecular mechanisms caused by tubulin structural polymorphism.

Starting from 1972, colchicine is known as the most useful drug for prevention of familial Mediterranean fever attacks. However, some patients do not respond to colchicine treatment, even taken in high doses. Despite the fact, that different hypotheses have been proposed, the molecular mechanisms of colchicine resistance are not completely clear. It is generally known, that colchicine binds ß-tubulin and inhibits microtubules polymerization. The ß-tubulin gene has SNPs, which lead to amino acid substitutions, and some of them are located in colchicine binding site (CBS). We have assumed, that this SNPs can affect tubulin-colchicine interaction and might be the reason for colchicine resistance. With this in mind, we modeled 7 amino acid substitutions in CBS, performed molecular dynamics simulations of tubulin-colchicine complex and calculated binding energies for every amino acid substitution. Thus, our study shows, that two amino acid substitutions in the ß-tubulin, namely A248T and M257V, reduce binding energy for approximately 2-fold. Based on this, we assume, that these amino acid substitutions could be the reason for colchicine resistance. Thus, our study gives a new insight into colchicine resistance mechanism and provides information for designing colchicine alternatives, that could be effective for colchicine resistant patients.

Inhibition of African swine fever virus infection by genkwanin.

African swine fever virus (ASFV) is the causative agent of an economically important disease of pigs for which no effective vaccines or antiviral drugs are available. Recent outbreaks in EU countries and China have highlighted the critical role of antiviral research in combating this disease. We have previously shown that apigenin, a naturally occurring plant flavone, possesses significant anti-ASFV activity. However, apigenin is practically insoluble in highly polar solvents and it occurs typically in derivative forms in plants. Here we screened several commercially available apigenin derivatives for their ability to inhibit ASFV Ba71V strain in Vero cells. Among them, genkwanin showed significant inhibition of ASFV, reducing viral titer from 6.5 ±â€¯0.1 to 4.75 ±â€¯0.25 log TCID/ml in a dose-dependent manner (IC50 = 2.9 µM and SI = 205.2). Genkwanin reduced the levels of ASFV early and late proteins, as well as viral DNA synthesis. Our further experiments indicated that genkwanin is able to inhibit ASFV infection at entry and egress stages. Finally, genkwanin displayed potent antiviral activity against highly virulent ASFV isolate currently circulating in Europe and China, emphasizing its value as candidate for antiviral drug development.

Genistein inhibits African swine fever virus replication in vitro by disrupting viral DNA synthesis.

African swine fever virus (ASFV) is the causal agent of a highly-contagious and fatal disease of domestic pigs, leading to serious socio-economic consequences in affected countries. Once, neither an anti-viral drug nor an effective vaccines are available, studies on new anti-ASFV molecules are urgently need. Recently, it has been shown that ASFV type II topoisomerase (ASFV-topo II) is inhibited by several fluoroquinolones (bacterial DNA topoisomerase inhibitors), raising the idea that this viral enzyme can be a potential target for drug development against ASFV. Here, we report that genistein hampers ASFV infection at non-cytotoxic concentrations in Vero cells and porcine macrophages. Interestingly, the antiviral activity of this isoflavone, previously described as a topo II poison in eukaryotes, is maximal when it is added to cells at middle-phase of infection (8 hpi), disrupting viral DNA replication, blocking the transcription of late viral genes as well as the synthesis of late viral proteins, reducing viral progeny. Further, the single cell electrophoresis analysis revealed the presence of fragmented ASFV genomes in cells exposed to genistein, suggesting that this molecule also acts as an ASFV-topo II poison and not as a reversible inhibitor. No antiviral effects were detected when genistein was added before or at entry phase of ASFV infection. Molecular docking studies demonstrated that genistein may interact with four residues of the ATP-binding site of ASFV-topo II (Asn-144, Val-146, Gly-147 and Leu-148), showing more binding affinity (-4.62 kcal/mol) than ATP4- (-3.02 kcal/mol), emphasizing the idea that this viral enzyme has an essential role during viral genome replication and can be a good target for drug development against ASFV.

Complex formation dynamics of native and mutated pyrin's B30.2 domain with caspase-1.

Pyrin protein is the product of the MEFV gene, mutations in which cause manifestation of familial Mediterranean fever (FMF). Functions of pyrin are not completely clear. The secondary structure of the pyrin is represented with four domains and two motifs. Mutations p.M680I, p.M694V, p.M694I, p.K695R, p.V726A, and p.A744S, which are located in the B30.2 domain of pyrin protein, are responsible for manifestation of the most common and severe forms of FMF. All the domains and the motifs of pyrin, are directly or indirectly, involved in the protein-protein interaction with proteins of apoptosis and regulate the cascade of inflammatory reactions, which is impaired due to pyrin mutations. It is well known, that malfunction of the pyrin-caspase-1 complex is the main reason of inflammation during FMF. Complete tertiary structure of pyrin and the effects of mutations in it are experimentally not studied yet. The aim of this study was to identify possible effects of the abovementioned mutations in the B30.2 domain tertiary structure and to determine their potential consequences in formation of the B30.2-caspase-1 complex. Using in silico methods, it was found, that these mutations led to structural rearrangements in B30.2 domain tertiary structure, causing shifts of binding sites and altering the interaction energy between B30.2 and caspase-1.

Molecular dynamics study of interaction and substrate channeling between neuron-specific enolase and B-type phosphoglycerate mutase.

Phosphoglycerate mutase (PGM) and enolase are consecutive enzymes in the glycolytic pathway. We used molecular dynamics simulation to examine the interaction of human B-type PGM (dPGM-B) and neuron-specific enolase (NSE). Specifically, we studied the interactions of 31 orientations of these enzymes by means of the effective energy function implicit solvation method. Interactions between active regions of the enzymes occurred preferentially, although the strongest interactions appeared to be between the back side of NSE and the active regions of dPGM-B. Cleavage of 2PG from dPGM-B was investigated, and the Ser(14)-Leu(30) loop of dPGM-B is suggested as a cleavage site and, likely, another entrance site of a ligand. Substrate channeling between the enzymes was observed when NSE with its active regions Leu(11)-Asn(16), Arg(49)-Lys(59), and Gly(155)-Ala(158) covered the Ser(14)-Leu(30) loop of dPGM-B. Analyses of the results make us believe that the channeling between PGM and enolase "benefits" from weak interaction. The probability of formation of channeling favorable complex is estimated to be up to 5%, while functional interaction between NSE and dPGM-B might be as high as 20%. NSE and dPGM-B functional interaction seems not to be isotype specific.